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Abstract 

We study a planar polynomial differential system, given by x = P(x,y), 
y = Q(x,y). We consider a function I(x,y) = ex.p{fi2(x)Ai(x, y)/Ao(x, y)} 
hi(x)Y\i = i(y — 9i{x)) a% , where gi(x) are algebraic functions of x, A\(x,y) = 
Hk=i(y - a k(x)), A (x,y) = l\j=i(y - 9j{x)) with a k {x) and gj(x) algebraic 
functions, Ao(x,y) and A\(x,y) do not share any common factor, h%(x) is a 
rational function, h(x) and h\(x) are functions of x with a rational logarithmic 
derivative and on £ C. We show that if I(x, y) is a first integral or an 
integrating factor, then I(x, y) is a Darboux function. A Darboux function is 
a function of the form f^ 1 . . . f p p exp{/i//o}, where fi and h are polynomials 
in C[x, y] and the Aj's are complex numbers. In order to prove this result, we 
show that if g(x) is an algebraic particular solution, that is, if there exists an 
irreducible polynomial f(x,y) such that f(x,g(x)) = 0, then f(x,y) = is 
an invariant algebraic curve of the system. In relation with this fact, we give 
some characteristics related to particular solutions and functions of the form 
I(x, y) such as the structure of their cofactor. 

Moreover, we consider Ao(x,y), Ai(x,y) and h,2(x) as before and a func- 
tion of the form := exp{/i2(x) A\{x, y)/Ao(x, y)}. We show that if 
the derivative of <3?(x, y) with respect to the flow is well defined over {(x, y) : 
Ao(x,y) = 0} then $(x,y) gives rise to an exponential factor. This exponen- 
tial factor has the form exp{R(x, y)} where R = h,2A\/Ao + B\/Bq and with 
Bi/Bq a function of the same form as /z^Ai/A). Hence, ex.p{R(x, y)} factor- 
izes as the product $(x, y) ^(x, y), for y) := exp{B\(x, y)/Bo(x, y)}. 
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1 Introduction 



In this work we consider planar polynomial differential systems as: 

x = P(x,y), y = Q(x,y), (1) 

where P(x,y) and Q(x,y) belong to the ring of real polynomials in two variables, 
M[x, y]. We will always assume that P(x,y) and Q(x,y) are coprime polynomials. 
We denote by d = max{degP, degQ} and we say that d is the degree of system (JTJ). 
Equivalently to system ([T}, we may consider the ordinary differential equation: 

dy_ = Q(x,y) (2) 
dx P{x,y)' 1 ' 

In order to shorten some formulae, we introduce the operator X associated to (Q): 

This paper is related with the study of the properties of a certain type of par- 
ticular solutions of system (JTJ). Given a function g(x), we say that y — g(x) is a 
particular solution of equation © if 

_ Q(x,g(x)) 
9[X) ~ P(x,g(x)Y (3) 

where g'(x) = dg(x)/dx. In particular, we are concerned with algebraic functions. 
We say that a function g(x) is algebraic if there exists a polynomial f(x, y) such that 
f(x,g(x)) = 0. As we prove in the next section, by using some algebraic results 
stated for instance in [H [TUl HZj, this polynomial f(x,y) can always be chosen 
irreducible and it is unique modulus multiplication by constants. 

The solutions of system may also be given in an implicit way. An invariant 
algebraic curve of system (JTJ) is an algebraic curve f(x,y) = satisfying 



df df 

P(x, y)^-(x, y) + Q(x, y)^-( x , y) 

ox ay 



= 0. (4) 

f(x,y)=o 



By using Hilbert's Nullstellensatz, cf. it can be shown that f(x,y) = is an 
invariant algebraic curve of system if, and only if, there exists a polynomial 
k(x,y) G R[x, y] satisfying: 

df df 

p(x, y)-Q^( x ' y) + Q( x i y)-Q-( x i y) = k ( x > y)f( x ^ v)- 

This polynomial k(x,y) is called the cofactor of the curve given by f(x,y) = 
and it can be shown that its degree is lower than or equal to d — 1. Invariant 
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algebraic curves, also denoted as Darboux polynomials by some authors, have been 
widely studied for their relation with integrability and some qualitative properties 
of polynomial differential systems, see for instance ^21 EH] • As a generalization of 
the notion of invariant algebraic curve, we can define an exponential factor. This 
concept has been firstly introduced by Christopher [3] and it is related with the 
notion of multiplicity of an invariant algebraic curve of system An exponential 
factor is a function of the form exp{h/ f } where h(x, y) and fo(x, y) G C[x, y], and 
X (exp{/i//o}) = k(x,y) exp{/i// } with k(x,y) a polynomial of degree at most 
d — 1. The following lemma is given in [3| and characterizes exponential factors. 

Lemma 1 The function exp{h/ f } is an exponential factor of system (JTJ) with co- 
factor k(x,y) if, and only if fo(x,y) = is an invariant algebraic curve of system 
(HJ) with cofactor k (x , y) and X{h) = k h + A; Jo- 
Invariant algebraic curves characterize the existence of first integrals for system ((TJ) 
belonging to a certain functional class. In order to properly state the known results 
about integrability using invariant algebraic curves, we need to consider complex 
algebraic curves f(x,y) = 0, where f(x,y) G C[x,y]. Since system JTJ is defined by 
real polynomials, if f(x,y) = is an invariant algebraic curve with cofactor k(x,y), 
then its conjugate f(x,y) = is also an invariant algebraic curve with cofactor 
k(x,y). Hence, its product f(x,y)f(x,y) G M.[x, y] gives rise to a real invariant 
algebraic curve with a real cofactor k(x,y) + k(x,y). In the Darboux theory of 
integrability, quoted in the forthcoming paragraph, we need to consider invariant 
algebraic curves defined by polynomials in C[x, y] since they play an essential role in 
the theory of integrability. We also notice that in M 2 the curve given by f(x, y) = 0, 
even if f(x, y) G M.[x, y], may only contain a finite number of isolated singular points 
or be the null set. 

A function of the form f^ 1 . . . fp p exp{h/ f }, where /j and h are polynomials in 
C[x, y] and the Aj G C, is called a Darboux function, see for instance [T2J[T3]. System 
(Q) is called Darboux integrable if it has a first integral or an integrating factor which 
is a Darboux function (for a definition of first integral and of integrating factor, 
see [2 ). The following lemma gives the relation between Darboux functions and 
invariant algebraic curves of a system (0) . 

Lemma 2 We consider a Darboux function T)(x, y) := f^ 1 . . . fp p exp{h/ f } such 
that 

dV dV 

P{x, y)-z-(x, y) + Q(x, y)—{x, y) = k(x, y)V(x, y), 
ox oy 

where k(x,y) is a polynomial of degree at most d — 1. Then, each fi(x,y) = 0, 
i = 1,2, ...,p is an invariant algebraic curve of system and exp{h/f } is an 
exponential factor of system (JTJ) . 
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The proof of this lemma is analogous to the proofs of Lemma 3 and Proposition 4 
in 0. 

We recall that V(x, y) is an inverse integrating factor of system (|TJ) if it is a 
function of class C 1 in some open set U of M. 2 and satisfies the following partial 
differential equation: 



We note that the function 1/V(x,y) is an integrating factor for system (JJJ in 14. 
The following result, which is a summary of some well known results, relates the 
existence of an inverse integrating factor in a certain functional class and the exis- 
tence of a first integral in another (possibly larger) functional class. The definitions 
of elementary function and Liouvillian function can be found in [T^]. Proposition 
E] shows that when considering the integrability problem we are also addressed to 
study whether an inverse integrating factor belongs to a certain given class of func- 
tions. As many authors have noted, see for instance [2], inverse integrating factors 
play a fundamental role in the integrability problem, not only because they charac- 
terize the functional class of a first integral but also because they usually belong to 
an easier functional class. For instance, quadratic systems of the form (|TJ) with a 
center at the origin always have an inverse integrating factor which is a polynomial 
of degree 3 or 5. Hence, we have that the characterization of centers for quadratic 
systems can be done by means of a polynomial instead of a first integral, which will 
be of Darboux type in a general case. 

Proposition 3 The following three statements hold. 

(a) System (JJJ has a Darboux first integral if, and only if, it has a rational inverse 
integrating factor. 

(b) // system (JJJ has an elementary first integral, then it has an inverse integrating 
factor of the form V(x,y) = (A(x, y)/B(x, y)) l ^ N , where N e Z and A, B G 



(c) System (JJJ) has a Liouvillian first integral if, and only if, it has a Darboux 
inverse integrating factor. 

The first statement of this proposition is proved in j2] and its reciprocal is proved 
in jUHS]. Statement (b) is proved in [T3] and the last statement is proved in [To] . 

It is clear, as shown in Lemma 2 of 0, that given an invariant algebraic curve 
f(x, y) = of system (JJJ), all the algebraic functions defined by it in an implicit way, 




C[x,y). 
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that is, all the functions g(x) satisfying f(x,g(x)) = 0, are particular solutions of 
equation (j2J). In this paper, among other results, the converse result is established, 
that is, given an algebraic particular solution y = g(x) of equation (j2J), we show 
that the irreducible polynomial f(x,y) such that f(x,g(x)) = gives rise to an in- 
variant algebraic curve f(x,y) = of system (JTJ). This fact is stated and proved in 
Theorem El of Section El We have noticed that a Darboux function may also contain 
exponential factors and this fact is necessary so as to characterize the Liouvillian 
integrability. Hence, exponential factors appear in a natural way when considering 
invariant algebraic curves. In this paper we consider algebraic particular solutions 
y — g(x) which come naturally from invariant algebraic curves and this relationship 
allows us to give an analogous to exponential factors but for algebraic particular 
solutions. This analogy is motivated and made clear in Subsection 13.21 In relation 
with this fact, we give some characteristics related to particular solutions such as 
the structure of their cofactor, which is given in Subsection 13.31 

In [6] and jH] an algorithmic method to determine, for system the possible 
existence of first integrals or integrating factors of the form I(x, y) = h(x) Y[ i=x {y — 
gi(x)) a * or I(x,y) = exp{h 2 (x) l\ r k=1 (y - a k (x))/ ]\ s j=1 (y - g s (x))} h(x) YlLiiv ~ 
gi(x)) ai , where gi(x) and (jj(x) are unknown particular solutions of equation (J2J), 
Oj 6 C are unknown constants, and dk(x), h(x), hi(x) and h 2 (x) are unknown 
functions, is given. In both cases, if all the particular solutions g-i{x) and (jj(x) are 
determined, which is expressed by the non-existence of a nonlinear superposition 
principle as described in (HI IE] , they are algebraic functions (see Proposition 7 in [6 j 
and Theorem 2 in [Hj). The algorithm, in this case, gives an alternative method to 
determine such type of solutions. In the case where all the gi(x), a k (x) and (jj(x) 
are algebraic functions, h(x) and h\{x) have a rational logarithmic derivative and 
h 2 (x) is a rational function, we show (cf. Propositions [TBI and IT71 Subsection I3.4|) 
that I{x,y) is a Darboux function. This result is one of the main goals obtained 
in the present work. Hence, using the algorithm described in [SJ |Hj , all the systems 
with a Liouvillian first integral can be found, as well as an explicit expression of a 
non-Liouvillian first integral when there is a nonlinear superposition principle. The 
present work is born as a complement to the results described in |3 IHJ - In these two 
works the integrability problem is studied in the particular case of the existence of 
a first integral of the described form I(x,y). As we have already stated, in [HI IE] 
it is shown that if when applying the described algorithm the function I(x, y) is 
completely determined then it only involves algebraic functions. In that case, in 
all the examples studied in [HIE], the function I(x,y) was a Darboux function. In 
fact, this is the general case: if when applying the algorithm described in jHIEj the 
function I{x,y) is completely determined then it is a Darboux function. The proof 
of this assertion is one of the objectives that we have achieved in the present work. 



5 



2 Some preliminary results on algebraic functions 
and polynomials 

In this section we give a summary of well known algebraic results which are needed 
in this paper. The definitions and proofs can be found, for instance, in the books 

nun ei . 

We always consider polynomials in M[x, y] which is a ring with the usual addi- 
tion and product operations of polynomials. Equivalently, we may consider the ring 
C[x, y]. One of the most important equivalence relations which can be defined in 
the ring M[x, y] is the divisibility relation. We say that the polynomial f\ G M.[x, y] 
divides f 2 £ IR.[or, 2/] , and we write fi \ f 2 , if there exists a polynomial k G M[x, y] 
such that f 2 = k f\. It can be shown that Wt[x, y] is a unique factorization domain 
(UFD). We recall that the unit elements of the ring M.[x, y] are the constants differ- 
ent from zero, i.e., M— {0}. We say that two elements /1, f 2 G M[x, y] are associates 
if there exists a unit element e such that f\ = e f 2 . An irreducible polynomial is a 
non-constant element f(x,y) G M.[x, y) which is only divided by its associates. 

If /i)/2 £ ^[^y] are such that for any (a, b) G {(x, y) G C 2 : f\(x,y) = 0} 
we have that f 2 (a,b) = 0, we will write f 2 (x, y)\f 1 ( x y ) =0 — 0- Given a polynomial 
/ G M[x, y], it defines a curve, denoted by f(x, y) = 0, which is the set {(x, y) G C 2 : 
f(x, y) = 0}. We say that p G C 2 is a point of intersection of two curves fi(x, y) = 
and f 2 (x, y) — if fi{p) = f2(p) = 0. An intersection point must be always counted 
as times as its multiplicity. The multiplicity of a point of intersection is a rather 
complicated notion, which can be found, for instance, in page 60 of ^7j. Intuitively, 
the multiplicity of a point of intersection takes into account the number of tangents 
shared by the two curves at that point. 

Bezout's theorem takes into account the degree of two curves and their points 
of intersection and relates them with the fact of having a common factor. 

Theorem 4 [Bezout] If fi(x,y) = and f 2 (x,y) = are two curves of degrees 
m and n, respectively, with more than mn intersection points, then there is an 
irreducible polynomial r(x,y) which divides both f\(x,y) and f 2 (x,y). 

We are going to recall some results on fractionary power series that can be 
found in P]. If x is any free variable over C, we denote by C((x)) the field of 
fractions of C[[x]], where C[[x]] is the ring of entire formal power series in x with 
complex coefficients. We recall that given a ring, which needs to be an integral 
domain, we can define its field of fractions as the smallest field containing it. The 
field of fractions is therefore obtained from the integral domain by adding the least 
needed to make of it a field, that is, the possibility of dividing by any nonzero 
element. Given n G N, we consider entire fractionary series s = ^ i>r aiX^ n where 
r G Z, <2j G C and min{i : <2j 7^ 0} > 0. An element s G C((x 1//n )) is of the form 
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s = X]i>r a i x% ^ n where r G Z, a, G C and the minji : a« 7^ 0} can be lower than zero. 
The elements of the ring Cp 1 /™]] are the entire fractionary series such that min{z : 
a>i 7^ 0} > 0. It can be shown (see pages 17 and 18 of [T]), that given a fractionary- 
series s = X]j>r a i x £ C((x 1 / n )) we can always take an equivalent series such that 
n and gcd{i : aj 7^ 0} have no common factor. In this case, we say that n is the 
polydromy order of the fractionary series s and we denote it by v{s). Let M(x,y) 
be a polynomial in y of degree N G N whose coefficients in ?/ are fractionary series 
in C((a; 1 / ni )), that is, we expand M(x,y) in powers of y: M{x,y) = J2iL Si(x)y l 
and we have that Si(x) G C((x 1//ni )), for i — 0, 1, 2, . . . , N. The polydromy order of 
M(x,y) is defined as the least common multiple of the polydromy orders of Si(x), 
i = 0, 1, 2, . . . , N, that is, u(M) = lcm{z/(s;) : i = 0, 1, 2, . . . , N}. 

The fractionary series s is said to be convergent if Xli>r a $ ^ as non-zero con- 
vergence radius, where t = x x l n . 

The ring C{x,y} is the ring of convergent power series in two variables and 
complex coefficients. The following result clarifies the structure of an algebraic 
function and it is stated and proved in [T] (page 26). 

Theorem 5 If f(x,y) G C[x, y], then there is a unique decomposition of the form 
f = ux r Y\i = i(y — 9i(%)), where r G N U {0}, u is an invertible power series of 
C{x, y} and gi(x) are convergent fractionary series. 

In the book pQ this result is stated for formal power series f(x,y) G C[[x, y]\. In 
the Corollary 1.5.6 of page 26 in pQ, it is stated that i is the height of the Newton 
polygon associated to f(x, y). Since we are considering a polynomial f(x, y) instead 
of a formal series, we deduce that £ is the highest degree in y of f(x, y) and, therefore, 
u is an invertible power series of C{x}. We define h(x) = ux r and we have that 
given f(x, y) G C[x, y] of degree £ in y, there is a unique decomposition of the form: 

e 

f = h(x) Y[(y- gi (x)), (5) 

i=l 

where h(x) G C{x} and gi(x) are fractionary series. Theorem 1.7.2 (page 31) of [1 
ensures that all the y-roots of f(x,y) are convergent. Hence, we deduce that an 
algebraic function g(x) is a convergent fractionary series. 

By definition, given an algebraic function g(x), we have that there exists a 
polynomial such that f(x,g(x)) = 0. If f(x,y) is not irreducible, then g(x) is a 
fractionary series appearing in the decomposition (0) of at least one of the irreducible 
factors of f(x,y). Hence, without loss of generality, we may always assume that 
f(x,y) is irreducible. Moreover, given g(x), this irreducible polynomial f(x,y) is 
unique (modulus multiplication by constants). This statement is clear from the fact 
that if fi(x, y) and f2^x, y) are two irreducible polynomials such that fi(x, g(x)) = 0, 
i — 1,2, then these two polynomials have an infinite number of points of intersection 
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(because g{x) is a convergent fractionary series) and, by Bezout's Theorem we 
have that f\{x,y) and f2{x,y) must be associates. 



3 The Main Results 

3.1 Particular algebraic solutions 

Theorem 6 Let g(x) be an algebraic particular solution of equation (J2J) and we 
call f(x,y) the irreducible polynomial satisfying f(x,g(x)) = 0. Then, the curve 
f(x,y) = is an invariant algebraic curve of system (JTJ). 

Proof. Let us denote by F(x,y) the polynomial in M.[x, y] defined by: 

df df 
F(x, y) := P(x, Vlg^ix, V) + Q{%, V)- 

We have that F(x,y) = and f(x,y) = intersect in all the points of the form 
(x, g(x)) by virtue of @. By Bezout's TheoremHJ we deduce that the polynomials 
f(x, y) and F(x, y) share a common factor, because they intersect in an infinite (con- 
tinuum) number of points. They intersect in all the points (x,g(x)) where the frac- 
tionary power series g(x) is convergent. Since f(x,y) is an irreducible polynomial, 
we have that f(x, y) divides F(x, y) in the ring of real polynomials. From this fact, 
we conclude that there exists a polynomial k(x, y) such that F(x, y) = k(x, y)f{x, y) 
and we get that f(x,y) = is an invariant algebraic curve of system | 

We recall the definition of invariant and of quasipolynomial cofactor stated in jH] . 

Definition 7 An invariant of ^ is a function <p(x, y) such that there exists a 
quasipolynomial cofactor M(x,y), where M(x,y) is a polynomial in one of the 
variables x or y of degree < m — 1 with m the degree of system (0) in that variable, 
satisfying P(d(j)/dx) + Q(d(f)/dy) = M<p. 

In case that the set of points in C 2 satisfying that <j)(x, y) — is not null, we have 
that (p is an invariant if, and only if, P(d(p/dx) + Q(d(f>/dy)\ ( f >=0 = 0. We say that 
= is an invariant curve in this case. 

These definitions are a generalization of the so called generalized cofactor in- 
troduced in |7j where a generalization of the Darboux integrability theory in order 
to find non-Liouvillian first integrals of system was presented. For the special 
invariant curve 4>(x,y) := y — g(x) = of (HJ), where g(x) is a particular solution 
of equation (j2J), a quasipolynomial cofactor always exists as it was established in 
[0]. Many examples of invariants with a quasipolynomial cofactor are given in |H] 
as well as a method to find first integrals, which are non-Liouvillian in general, for 
certain families of systems. 
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Proposition 8 [6 J A particular solution g(x) of equation (J2J) always has a unique 
associated quasipolynomial cof actor of the form M(x,y) = fc TO _i(x)y m ~ 1 + ••• + 
kx(x)y + k (x), where m is the degree in y of system ([!}. 

In case that g(x) is algebraic then each ki(x), i — 0, 1, 2, . . . , m — 1 is a rational 
function in x and g(x) with coefficients in C. 

The second part of this proposition is deduced from the proof of the first part as 
given in [Hj. 

In case g(x) is an algebraic function, by Theorem we have that it is a frac- 
tionary series. Since ki(x) is a rational function in x and g(x) with coefficients in 
C and g{x) e C((x 1 / n )), for a certain natural n, then its quasipolynomial cofactor 
M(x,y) G C((x^ n ))[y\, that is, each of the functions fcj(x), i = 0, 1, 2, . . . , m — 1 
is a fractionary series with a polydromy order that divides the polydromy order of 
g(x). Next lemma shows that in case M(x,y) is a polynomial, then g(x) must be a 
rational function. 

We say that an equation (J2J) is linear if it is of the form 



where m (x) and mi(x) are functions of x. If equation (j2J) does not take this form, 
we say that it is non-linear. 

Lemma 9 Let g(x) be a particular solution of a rational non-linear equation (J2J). 
If the quasipolynomial cofactor M(x,y) related to g(x) is a polynomial, then g(x) is 
a rational function. 

Proof. We have that M(x, y) = k m _i(x)y m ~ 1 + • • • + ki(x)y + k (x) is a polyno- 
mial, so, kj(x) is a polynomial in x for all j = 0, 1, 2, . . . , m — 1. Let us expand the 
polynomials P(x,y) and Q(x,y) in powers of y: P(x,y) = Po(x) + pi(x)y + . . . + 
Pm(x)y m , Q(x,y) = qo(x) + qi(x)y + . . . + q m (x)y m , where Pj(x) and qj(x) are the 
polynomials in x corresponding to the coefficients of degree j in y of P(x,y) and 
Q(x,y), respectively. Since g(x) is a particular solution of equation ©, we have 
that —P(x,y)g'(x) +Q(x,y) = M(x,y)(y — g(x)). Equating the coefficients of order 
j my, we deduce that: 

p (x)g'(x) - k (x)g(x) = q (x), (eg ) 

Pj(x)g'(x) - kj(x)g(x) = q^x) - k,-x(x), (eqj) j = 1, 2, . . . , m - 1, 

p m (x)g'(x) = q m (x) - k m -i{x) (eq m ). 

If k (x)pj(x) —p (x)kj(x) ^ for some j = 1, 2, . . . , m — 1, we can equate g(x) from 
the equations (ego) and (eg,) and we get 



dy_ 

dx 



mi(x) y + m (x), 



Po(x)qj(x) - qo(x)pj(x) - p Q {x)k j _ 1 {x) 
k (x)pj(x) -p {x)kj(x) 
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which is a rational function. If k Q (x)pj(x) — po(x)kj(x) = Oforallj = 1,2, ... ,m— 1, 
we deduce that (x) = fcj(a;)L 1 (a;) andpo^) — ^0(^)^1 (^) f° r all j = 1,2, ...,m— 1, 
where Li(x) is a rational function in x. 

If p m (x)A;o(a;) ^ 0, from the first and the last equations we deduce that 

, s _ Po(^)(gm(^) ~ fcm-i(^)) -Pm(x)g (a:) 

which is a rational function. 

We can also try to equate g(x) from the equations (eq m ) and (eg 3 ) for some 
j = 1, 2, . . . , m — 1. If p m (x)kj(x) ^ 0, then 

_ q m {x)pj(x) -p m {x)qj{x) + p m (x)kj^i(x) - k m ^i(x)pj(x) 

p m (x)kj(x) 

which is a rational function. 

In case that ko(x)pj(x) — po(x)kj(x) = 0, p m (x)ko(x) = and p m (x)kj(x) = 
for all j = 1, 2, . . . , m — 1, first assume that M(x, y) = 0. Then y — g(x) would be 
a first integral for system (JTJ), which means that Q(x, y) — P(x, y)g\x) = 0. Hence, 
Q(x, y) and P(x, y) are such that equation © is a linear one, in contradiction with 
our hypothesis. So, we conclude that p m (x) = 0. We also have that P(x, y) = 
Li(x)M(x, y), that is, Pj(x) = fcj(x)L 1 (a;) for j — 0, 1, 2, . . . , m — 1. From equation 
(eq m ) we have that q m (x) = k m _i(x) and the other equations read for: 

k (x)(L 1 (x)g'(x) - g{x)) = q (x), (eq ), 

k j (x)(L 1 (x)g'(x) - g{x)) = qj(x) - k^ x {x), (egj-), 

for j = 1,2, . . . ,m — 1. Equating the factor (Li(x)g'(x) — g(x)) from equation 
(eq' ) and (eq'j) we deduce that kj(x)q (x) — k (x)qj(x) + k (x)kj_i(x) = for all 
j = 1, 2, . . . ,m — 1. We write (fo^) — -^0(^)^0(^)5 where I/o(^) is a rational function, 
and we have that k (x)(kj(x)L (x) —qj(x) +kj_i(x)) = for all j — 1, 2, . . . , m— 1. 
If k (x) = 0, then P(x, y) and Q(x, y) share the common factor y. If ko(x) ^ 0, then 
we have that qj{x) = kj(x)L (x) + kj-\(x), from which we deduce that Q(x,y) = 
{y + L {x))M{x,y). 

Unless M(x,y) is a real number (different from zero), we have that P(x,y) and 
Q(x,y) share a common factor, in contradiction with our hypothesis. Therefore, 
after a rescaling of the time, the only systems of the form (0) with a particular non 
rational solution g(x) with a polynomial cofactor (in fact, the cofactor is the real 
number 1) are of the form : x — L\(x), y — y + -^0(^)5 which give rise to a linear 
equation. | 

In the same way as in Lemma |H1 we have that if g(x) is an algebraic particular 
solution of a nonlinear equation (J2J) with polydromy order v(g), then its quasipoly- 
nomial cofactor must have the same polydromy order. 
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Lemma 10 Let g(x) be an algebraic particular solution of a non-linear equation 
(J2J). Its quasipolynomial cofactor M(x,y) has the same polydromy order as g(x). 

Proof. We have that if g(x) has polydromy order v(g), then the polydromy order 
of its quasipolynomial cofactor M(x,y), z^(M), divides it, that is, z/(M)|z/(g). This 
is because M(x,y) is a polynomial in y and a rational function in x and g(x) as it 
has been stated in Proposition |H| 

The fact that both polydromy orders coincide is a corollary of the previous proof 
of Lemma EJ The reasonings are the same since we can equate g(x) in terms of the 
coefficients of M(x, y) unless the equation is of linear type. | 

REMARK 11 If we have a linear equation of the form dy/dx = rrii(x) y + mo(i), 
we may have a nonrational particular solution with a polynomial cofactor. For 
instance, taking m (x) = and m>i(x) = 1, we have that g(x) := e x is a particular 
solution with cofactor 1. In the same way, the linear equation dy/dx = 3y/(2x) 
has the algebraic solution y = x 3 ^ 2 , whose polydromy order is 2, with the polynomial 
cofactor 3, whose polydromy order is 1. 

3.2 On the invariants giving rise to exponential factors 

In the Darboux theory of integrability not only invariant algebraic curves are con- 
sidered, but also exponential factors, as we have stated in the introduction. Expo- 
nential factors appear in a natural way from the coalescence of invariant algebraic 
curves, as it is explained in j^j. This statement means that if we have a polyno- 
mial system (0) with an exponential factor of the form exp{h/f }, then there is 
a 1 -parameter perturbation of system (JT|), given by a small e, with two invariant 
algebraic curves, namely fo = and fo + eh = 0. Hence, when e = 0, these two 
curves coalesce giving the exponential factor exp{h/fo} for the system with e — 0, 
as well as the invariant algebraic curve /o = which does not disappear. 

In this context, and taking algebraic particular solutions into account, the fol- 
lowing question arises: which kind of function appears with the coalescence of two 
algebraic particular solutions? In the next Proposition we show that the natural 
generalization of algebraic particular solutions in this framework is a function of the 
form: $(x,y) := exp{h 2 (x) A^x, y)/A (x, y)} where At(x,y) = YYk=i(v ~ a k{ x )) 
and A Q (x, y) = YYj=i(y ~ 9j( x ))i with a^^x) and (jj(x) algebraic functions, k = 
1,2, ... , r and j — 1, 2, . . . , s, Ai(x, y) and A (x, y) do not share any common fac- 
tor, and h%{x) is a rational function in x. For convention, if r = or s = 0, we 
mean that Ai(x,y) or A (x,y) takes a constant value, respectively, which we may 
assume to be equal to 1. 

Since by Proposition |H1 we always have that a particular solution given by 
4>(x,y) = y — g(x) has a quasipolynomial cofactor, the property of being an in- 
variant for is given by associating to it a quasipolynomial cofactor M(x,y) 
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which is a polynomial in y of degree at most m — 1, where m is the degree in y of sys- 
tem Moreover, we will assume that M(x,y) is well defined over A (x,y) = 0, 
that is, M(x,Qj(x)) is a real function of x for all j = 1,2, ...,s. The following 
proposition gives the characterization of this fact. 

Proposition 12 We consider a function Q(x,y) := exp{h 2 (x) A\(x, y)/A (x, y)} 
where A x (x,y) = Yil=i(y ~ a k{x)) and A (x,y) = Y\ s j=1 {y - 9j(x)), with a k (x) and 
(jj(x) algebraic functions, k = 1,2, ...,r and j = 1,2, ...,s Ai(x,y) and A (x,y) 
do not share any common factor, and h 2 (x) is a rational function in x. Assume 
that there exists a function M(x,y), which is a polynomial in y of degree at most 
m — 1, where m is the degree in y of system such that: 

P(x, V)^^- + vf-^f 1 - = M(x, y) *(x, y). (6) 

We assume that M(x, gj(x)) is a real function of x for all j = 1,2, ...,s. We 
denote by ki(x) the coefficient of degree i iny of M(x, y), that is, M(x, y) = k (x) + 
k\(x)y + k 2 (x)y 2 + ... + k m _i(x)y m ' 1 . Then, 

(i) Each one of the algebraic functions gj(x), j = 1,2, ... ,s is a particular so- 
lution of equation (J2J). We denote by Mj(x,y) its associated quasipolynomial 
cofactor. 

(ii) The following identity is satisfied: 

x(h 2 {x)A 1 (x,y)^ = (^Mjfay) ) h 2 (x)A 1 (x,y) + M(x,y) A (x,y). (7) 



(iii) Each one of the functions ki(x) is rational in x and rational in ak(x) and 
gj(x), with k = 1,2, ... ,r and j = 1, 2, . . . , s. 



Proof. We have that X {<&(x, y)) = M(x,y)§(x,y), from which we deduce the 
following identity: 

X (h 2 (x) A x {x, y)) A (x, y) - h 2 (x)A 1 (x, y) X (A (x, y)) = M(x, y) A (x, yf . (8) 

Since A\[x, y) and A Q (x,y) do not share any common factor and A (x, (jj{x)) = 
for j = l,2,...,s, then, from equation (JHJ), we have that X (A (x, y)) \ y =gj x ) 
0. We notice that here we are using that M(x,gj(x)) is a real function of x for 
all j = 1,2, ... ,s. Let us call A Q . (x, y) := ni=i i^j(y ~ 9i( x )) and we have that 
A (x,y) = A 0j (x,y)(y - gj(x)). Then, 

X (A (x, y)) = X (A 0j (x, y)) (y - ~g 3 {x)) + A 0j (x, y) (-~g](x)P(x, y) + Q(x, y)) . 

12 



From X (A (x,y)) \ y= ^.^ = 0, we deduce that — g'j(x)P(x, (jj(x)) + Q(x, cjj(x)) = 
and we conclude that (jj(x) is a particular solution of system (JTJ), as stated in (i). 

In order to prove (ii), we consider the quasipolynomial cofactor Mj(x,y) associ- 
ated to y — (jj(x), whose existence is ensured by Proposition |H1 Then, we have that 

X(A ) = (EJ=i M j) A o- Hence, identity © reads for: 

x(h 2 {x)A 1 (x,y) ] jAo{x,y)-h2{x)A 1 (x,y) ^^Mj(x,y)\ A (x,y) = 

= M{x,y)A {x,yf. 

This identity coincides with (J7|), after dividing both members by A (x, y). 

Finally, to prove (iii), we observe that equating the coefficients of the same 
degree in y from identity ((7|), we deduce that each ki(x) is a rational function of x, 
gj(x), ak(x) and a' k (x). Since each afc(x) is an algebraic function, we can consider 
f ak (x,y) as the irreducible polynomial in C[x, y] such that f ak (x,ak(x)) = 0. We 
can derive this last identity with respect to x and we deduce that: 

dfg k (x,a k (x)) dfg k (x } a k (x)) _ 
dx + dy ak[X) = °- 

Hence, we can substitute the value a' k (x) appearing in ki(x) by the rational expres- 
sion in x and afc(x): — [df ak (x, ak{x))/dx] / [df ak (x, ak(x))/dy]. Therefore, each 
ki(x) is a rational function of x, Qj(x) and a>k(x), as we wanted to prove. | 

As Theorem |U1 shows, an algebraic particular solution y — g(x) recovers the 
invariant algebraic curve to which it is related. That is, we have that f(x,y) = is 
an invariant algebraic curve of system f^Q) if, and only if, all its y-roots are algebraic 
particular solutions of equation (j2J). The fact of being particular solutions implies 
that each one of them has an associated quasipolynomial cofactor. We notice that 
all the y-roots appear in the factorization f(x,y) = h{x) Y\^ = i{y — 9i{x)) described 
by ©• We would like to have an analogous to this statement but for exponential 
factors and this is what is given in the next Theorem 1131 

In this context, if we have an exponential factor given by exp{R(x, y)}, where 
R(x,y) is a rational function, i.e. R(x,y) G WL(x,y), we would like to write a 
factorization for it analogous to ©• We notice that given any two functions 
Ri(x,y), R2^x,y) with the property that R\(x,y) + R2(x,y) = R(x,y), we have 
that exp{R(x,y)} = exp{Ri(x, y)} ■ exp{R 2 (x,y)}, which is a "factorization" of 
exp{R(x,y)}. However, not all these factorizations are useful to our purposes be- 
cause it can be shown that exp{Ri(x, y)} and exp{R 2 (x, y)} do not need to have 
an associated quasipolynomial cofactor. So, they are not invariants for system 
as the exponential factor exp{R(x,y)} is. The following example exhibits this fact. 
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Example. Let us consider the following cubic system with the invariant straight 
line y = : 

x = (2x + y)(l + x) + 2x 2 y + y 3 , y — y(l+x + xy). (9) 

This system has the invariant y) := eyLY>{\fx / 'y} with the quasipolynomial co- 
factor M(x,y) = (1 + x + y 2 )/{2^/x). As it will be proved in Theorem FEU the 
existence of this invariant implies the existence of the following exponential fac- 
tor: F(x,y) := exp{(x + y)/y 2 } with the cofactor k (x,y) = y — x. We notice 
that not all the functions R\,R% satisfying Ri + R 2 = (x + y)/y 2 give rise to an 
invariant. For instance, if we take R\ := x/y 2 and R 2 '■— 1/y, the following easy 
computation shows that exp{i?!} is not an invariant of system ©, and neither 
exp{i? 2 } is. We note that, taking X as the vector field associated to system (0), we 
have that Af(exp{i?!}) = (y + (1 + x)/y) exp{Ri} and since the rational function 
(y + (1 + x)/y) is not well-defined over y = 0, we deduce that exp{_Rx} is not an 
invariant. However, if we consider y) := exp{(x + y — yyfx)/y 2 }, we have that 
y) ■ i/) = ?/) and both y) and ?/) are invariants of system (jHJ). 

The following theorem shows that if y) := exp{/i 2 (^) ^4i(^, y)/Ao(x, y)}, 
as described in Proposition H21 is such that it has an associated quasipolynomial 
cofactor, then there exists another function ^/(x,y) := exp{i? 1 (x, y)/B (x, y)}, of 
the same type as &(x,y), such that the product y)^(x, y) gives rise to an 
exponential factor exp{R(x, y)} for system This sentence means that a factor- 
ization of exp{R(x, y)} giving invariants for system is obtained by the product 
y)^(x, y). Therefore, we are able to recover exp{R(x, y)} from one of its fac- 
tors Q(x,y) and, moreover, since exp{i?(x, y)} appears by the coalescence of two 
invariant algebraic curves, we conclude that needs to appear by the coa- 

lescence of algebraic particular solutions since it is formed by a product of them. 
This fact is the result that we were targeting to: we wanted to exhibit the anal- 
ogy between the generalization of invariant algebraic curves to exponential factors 
with the generalization of algebraic particular solutions to invariants of the form 
$(x, y) := exp{h 2 (x) A x (x, y)/A (x, y)}. 

Theorem 13 Assume that the function <&(x,y) := exp{h 2 (x) Ai(x,y)/Ao(x,y)} 
has a quasipolynomial cofactor M(x,y), that is, X($(x,y)) = M(x,y) $(x,y) and 
assume that M(x,y) is well defined over {(x,y) : A (x,y) = 0}. Then, there exist 
quasipolynomial functions Bo(x,y) and Bi(x,y), which are polynomials in y and 
algebraic in x, such that R := h 2 A\/Aq + Bx/B is a rational function in x and y 
and exp{R(x,y)} is an exponential factor of system 

Proof. We consider (Pa(x, y, e) := A Q (x, y) +eh 2 (x) Ai(x, y) which is a polynomial in 
y whose coefficients are algebraic functions. We can compute the sequence of powers 
4>a{x, y, s) J+ , for each natural number j and the coefficients of all these polynomials 
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in y are algebraic functions, combination of those of 4>A(x,y,e). Therefore, there 
exists a natural number N such that <Pa(x, V, s) N+1 is a linear combination of all 
the previous powers and the coefficients of this combination are powers of x, y and 
e. This fact is due to the fmiteness of the algebraic extensions given by the y- 
roots of 4>a(x, y, e). In this way, we have that there exists an irreducible polynomial 
V(x, y, e) in x, y and e, such that each of the y-roots of 0a (x, y, e) is an y-root of 
V(x,y,e). This polynomial V(x,y,e) is the minimal polynomial of the y-roots of 
4>A{x,y,e), and it can be always taken irreducible. We expand V(x,y,e) in powers 
of e and we denote by Ri(x, y) its coefficient of degree % in e, which is a polynomial 
in x and y. Let us consider the quotient V(x, y, s)/(f>A(x, y, e) which is a polynomial 
in y, denoted by 4>b(%, y, s)- We expand 4>b (%, y, e) in powers of e and we denote by 
Bi(x, y) its coefficient of degree i in e, which is a polynomial in y. Since V = 4>a4>b, 
we deduce that Ro = A B and Ri = A Bx + h 2 AxB , equating the coefficients of 
e° and e 1 . Thus, R := h,2Ai/A + B\j B = Rx/Ro is a rational function in x and 
y. In case R\ = all the reasoning works just taking the first Ri ^ with i > 0. 

We only need to see that exp{R(x, y)} is an exponential factor of system (JTJ). We 
have, from Proposition that A is the product of algebraic particular solutions 
y — (jj(x). Hence X (A (x,y)) = M (x,y)A (x,y), where M (x,y) is the quasipoly- 
nomial cofactor X^=i Mj(x, y), which is a polynomial in y of degree at most m — 1, 
m being the degree in y of system (JTJ). By the irreducibility of V(x,y,e), we de- 
duce that Rq(x, y) is a power of the lowest degree polynomial containing as y-roots 
all the (jj{x). Hence, by Theorem EJ we have that Ro(x,y) = is an invariant 
algebraic curve of system (JTJ), that is, X(R ) = k R , where k (x,y) is a poly- 
nomial in x and y. Moreover, from Rq = AqBq and X (Aq) = MqAq we deduce 
that X (Bq) = (ko — Mq)Bq. By Lemma ^ we only need to show that there is a 
polynomial k Q (x,y) such that X (i?i) = k Q Ri + k R . We recall that from Propo- 
sition^] we have that X (h 2 Ai) = M Q h 2 Ai + M A . We consider the polynomial 
G(x,y) := X (R x (x,y)) - k (x,y)Rt(x,y). We have: 

G = X(R 1 )-k R 1 

= X (A ) B 1 + A X (fix) + X (h 2 A x ) B + h 2 A 1 X (B ) - k (A B X + h 2 A x B Q ) 
= M A B 1 + A X (50 + (MohA! + MA )B + h 2 A 1 (k - M )B 

- k (A B 1 + h 2 AtB ) 
= [MqBx + X (B x ) +MB - k Bx] A . 

We deduce that G(x, Qj(x)) = for all the cjj(x) appearing in A Q . We note that here 
we are using the hypothesis that M(x, y) is well defined over {(x, y) : A (x, y) = 0}. 
Since R (x,y) is the lowest degree polynomial with this property and G(x,y) is a 
polynomial, we have, by Bezout's Theorem, that Ro(x,y) is a divisor of G(x,y) 
in the ring of polynomials C[x, y]. So, there exists a polynomial ko(x,y) such that 
G(x, y) = k Q (x, y)Ro(x, y) and, thus, X (Rx) = k Rx + koRo. Therefore, exp{Rx/ Ro} 
is an exponential factor of system | 
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3.3 On the structure of the quasipolynomial cofactor 

The following proposition gives the form of the quasipolynomial cofactor associated 
to an invariant of the form I(x,y) = h(x) Yii=i{v ~ 9i( x )) ai - 

We define as m the degree of system in the variable y and we expand the 
polynomials P(x,y) and Q(x,y) in this variable: 

m m 



Proposition 14 Let I(x,y) = h(x) Yli=i(y ~ 9i{ x )) ai be an invariant of system 
(JTJ) with an associated quasipolynomial cofactor k(x,y) := k (x) + ki(x)y + . . . + 
k m ^i(x)y m ~ 1 . Then, p m (x)h'{x) = and 

h'(x) m ( a' s _(x) \ 

k A x ) = Pj( x ) j^j + Yl y a s-u+i)( x ) Qs( x ) - (g~- j) Ps( x )j » ( 10 ) 

/or j = 0, 1, 2, . . . , m — 1, where 

t 

°*k x ) = Y^gu^), for K = 0, 1,2, . ..,771. 
y=l 

VFe notice that the first term in (jlOj) does not appear if h(x) is a constant or the 
last term does not appear if p m ( x ) is zero. 

Proof. As it has been shown in [S], each of the factors y — gi(x) involved in the 
expression I(x,y) is a particular solution of system (JTJ). Moreover, a particular 
solution y — g(x) has a related quasipolynomial cofactor of degree at most m — 1 in 
y, therefore we have that k(x, y) is a polynomial in y of degree at most m — 1. We 
will deduce each expression of k{(x) from the identity: 



dl dl 

P(x, y)^-( x , y) + Q( x , v)-k-( x , y) = k(x, y)I(x, y). 
ox ay 



We compute the partial derivatives of I(x, y) in (JTTJ) and we divide each member of 
the resulting expression by I(x,y). We obtain: 



P(x,y) 



h'(x) _ y> a v g' v {x) 
K x ) ^iV-9u{ x ) 



+ Q( x ,y) 



E 



\y- 9u( x ) 



k(x,y). 
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Then, we deduce: 



P{x,y) 



h'{x) 
h(x) 



+ 



(12) 



The expression fJ12j) is an identity of polynomials in y of degree m + i. Let us 
consider the equality of coefficients of degree m + 1 in y: p m (x)h'(x)/h(x) = 0, 
which implies that either p m (x) = or h(x) is a constant function. 

We define eg, as the equation resulting from identifying the coefficients of y l in 



both members of equation ([12" 



0, 1, 2, . . . , m + £ — 1. From eqe +m -i, we can 



equate the expression of k m -i(x). Once we know this expression, from eqe +m -2, we 
can equate the expression of k m ^ 2 (x). Once we know this function, from e^ +m _3 
we can equate the expression of /c m _ 3 (x), and so on. Hence, in a recursive way, from 
eqi + j, we equate kj(x), where j = m — 1, m — 2, m — 3, . . . , 2, 1, 0. It can be shown 
by induction that these expressions are given by (JTUJ). We notice that 



(TAX 



K 



22 



ayg. 



K-l, 



x)gl(x) 



for K = 0, 1, 2, . . . , m. 

We substitute the values of kj(x) given in (fTUj) in equation (fT2~j) and we deduce 
that a polynomial of degree at most £ — 1 in y must be zero. We denote by Pol(y) 
this polynomial in y, which is also a function of x, h(x) and gi(x), g^fa), ■ ■ ■ ,ge( x )- 
We fix an index i such that 1 < i < £ and when we substitute y by gi(x) in Pol(y), 
we get that: 



Pol(gi(x)) = Yl (9i(x) -9b(x)) Q(x,gi(x)) - P{x,g i {x))g' i {x) 

Since each gi{x) is a particular solution of system (JTJ), we deduce that each ^j(x) 
is a different root of the polynomial Pol(y). Then, we have £ different roots of a 
polynomial of degree at most £ — 1, so Pol(y) is the null polynomial. We conclude 
that once we have substituted the cofactor as defined by (fTUjl in equation (|T2|) we 
have that this equation is satisfied. | 

In the same way as in Proposition El we are going to give the form of the 
quasipolynomial cofactor associated to an invariant of the form 



I(x, y) = exp < h 2 (x 



M{x,y) 
A (x,y) 



(13) 
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where we have defined A t (x,y) := lT k=1 (y - a k (x)) and A (x,y) := Y[) =1 {y -9j{x)). 

The first thing that we notice is that, without loss of generality, we can assume 
that r — s. In case that r / s, we can consider the change of variable y = 1/z. 
After this change (and a reparameterization of the time variable if necessary), we 
obtain another polynomial system (0) and the function / is transformed to another 
one with the same structure but with r = s. This last assertion is clear from the 
following equality: 



f nj=i (i - ft-(^) n;=i(* 



where ft^O^) := ^2(^)( — l) s ~ r ni=i a fe( a; )/ / llj=i ^ we ensure the structure 

of a quasipolynomial cofactor for the system with variables (x,z), we deduce its 
structure for the system with variables (x, y) just undoing the change and the repa- 
rameterization of the time, if it has been done. So, without loss of generality, we 
can assume that r = s. In fact, an analogous proof can be done for the case r ^ s 
but as it involves many computations which are not far from the ones that we are 
showing, we have preferred to avoid the case r ^ s by means of the change of 
variable y — 1/z. 

Proposition 15 Let I(x, y) be an invariant of system (JIJ) of the form (pHIj) with 
r = s and with an associated quasipolynomial cofactor k(x,y) := k (x) + k\{x)y + 
. . . + k m _i (x)y m ~ 1 . Then, p m (z)[/i' 1 (x) + hi(x)h' 2 (x)] = an d 

= (|H + + E «(*) - jr^J Mx) ) ' (14) 

for j = 0, 1, 2, . . . , m — 1, where either the first term does not appear in case that 
h\(x) = cexp{— ti2(x)}, with c G C and c ^ 0, or the last term does not appear if 
p m (x) = 0. The functions a K (x) are defined in the following way. Given k G N we 
consider the set of indexes 



r 



J K := < (ei, e 2 , . . . , e r , h, i 2 , • • • , v) : E Ek + E ^ = K ' e i G {°' 1 }' ^ G N U W 

{ k=l j=l 

and we have that: 

t 

~a K : = £ a u9 : + Kh 2 ^_ 1 yi+e a +..-+e r+ i ^ . . . a e r ~^n ...~ g 

v=l J K 

for k = 0, 1, 2, . . . , m. 



7, % T 

s ' 
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Proof. We have that if I(x,y) is an invariant, then each one of the functions gi{x) 
and (jj{x) are particular solutions and the identity (J7J) is satisfied. We are going to 
deduce the expressions of kj(x) only assuming that the following identity is satisfied: 

dl dl 



We multiply (jl5j) by 



A (x,^exp{- ^l|^ }n fa - ftM )^ 



so as to get an identity of polynomials in y of degree m+£+2r, where m is the degree 
of system (JIJ) in y. The equality of the coefficients of highest degree in y, that is 
ym+e+2r^ gj ves p m (x)[h[(x) + hi(x)h' 2 (x)] = 0. This condition gives two possibilities: 
either p m (x) = or hi(x) = cexp{— h 2 (x)}, with c G C and c ^ 0. Since the proof 
for both cases is analogous, we are going to follow them simultaneously. It has 
been shown in (8 i that k(x,y) is a polynomial in y of degree at most m — 1 so we 
can collect it in this variable: k(x,y) = Yll=i ki{x)y l . The equality of coefficients 
of degree j + £ + 2r gives us the expression of kj(x) in a recursive way. We first 
compute fc m _i from the equality of coefficients of y m - 1 + £ + 2r ; once we have this one, 
we compute k m _ 2 from the equality of coefficients of y m ~ 2 + e + 2r anc l so on. Some 
tedious computations show that these expressions are the ones given in (JHJ). 

We substitute the given expressions of kj(x) in the equality (fTKj) and from the 
equality of the coefficients of y e+2r ~ l we can compute the function h,2(x). Now we 
have a polynomial in y of degree at most i + 2r — 1 in y that must be identically 
zero. We denote by Pol(y) this polynomial in y. In the same way as in the proof of 
Proposition HH we substitute the y variable by each one of the £ + 2r functions gi(x), 
Qj{x) and ak(x) and we deduce that Pol(gi(x)) = because gi{x) is a particular 
solution for i = 1, 2, . . . , £, Pol(gj(x)) = because gj(x) is a particular solution for 
j = 1, 2, . . . , r and Pol(cik(x)) = because relation ((7|) is satisfied for k = 1, 2, . . . , r. 
Therefore, we have £ + 2r different roots of a polynomial of degree at most £ + 2r — 1, 
so Pol(y) is the null polynomial. We conclude that once we have substituted the 
cofactor defined by ()14|) and the function h 2 {x) in equation (fTo^) we have that this 
equation is satisfied. | 

Example. We include an example of the form of the quasipolynomial cofactor 
of a function I(x,y) as given by Proposition Hoi deduced from the equation (JHJ). 
We consider the following system: 

x = y + y + x + Ayx , y = —x — 2x + 2xy , (16) 

which has the exponential factor I(x,y) = exp{(2y — l)/(x 2 + y 2 )} with cofactor 
k(x,y) := —Ax. We are going to deduce the value of this cofactor by using the 
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formulas stated in (f*Hj) of Proposition [T^l We notice that this function I(x,y) has 
r = 1 < s = 2, so we need to perform the change of variables y = 1/z. The 
transformed system needs a reparameterization of the time, meaning multiplying it 
by z 2 , in order to have a polynomial system. Then, the transformed system reads 
for: 

x = 1 + (1 + Ax 2 )z + x 2 z 2 , z = -2xz 2 + (x + 2x i )z i , (17) 

which has the exponential factor I(x, z) = exp{z(2 — z)/(x 2 z 2 + 1)} with cofactor 
k(x, z) = —Axz 2 . This cofactor coincides with the transformation of the cofactor in 
coordinates (x, y) after the reparameterization of time. Using the notation described 
in PropositionlTolwe have that, for system (|17|) . m = 4, r = s = 2, i = 0, a\(x) := 0, 
a,2(x) := 2, h 2 (x) := — 1/x 2 , gi(x) := —i/x, g 2 {x) := i/x and h\(x) := 1, where 
i = y— 1. We recall that pt(x) corresponds to the coefficient of degree % in z of the 
polynomial defining x and qi{x) corresponds to the coefficient of degree % in z of the 
polynomial defining z. We have that p<±(x) = 0, so the first assertion of Proposition 
Hoi is satisfied. We compute h' 1 (x)/hi(x) +h' 2 (x) which gives 2/x 3 and the values of 
the a K (x) are equal to: &o(x) = 0, ai(x) = —2/x 2 , 5"2(x) = —2/x 4 , 03(2) = 6/x 4 , 
(ji(x) = 4/x 6 . An easy computation shows that the formulas written in (jl4j) give 
ko(x) = ki(x) = ks(x) = and k 2 (x) = — Ax, which corresponds to the given value 
of the cofactor k(x,z) = —Axz 2 . 



3.4 Darboux functions obtained from invariants 

We consider a Darboux function I(x, y) which is the product of invariant algebraic 
curves up to complex numbers. Then, by computing |/-roots, it can be expressed 
in the form I(x,y) = h(x) Yl i==1 (y — 9i{x)) ai , where G C, gi{x) are algebraic 
particular solutions and h(x) is such that its logarithmic derivative is a rational 
function. We recall that the logarithmic derivative of a function h(x) is the quotient 
h'{x)/h{x). 

In this subsection, we are concerned with the converse of this problem, that is, 
we give the conditions that a function of the form I(x, y) = h(x) Yl i= i(y — gi{x)) ai 
must satisfy in order to be a Darboux function. 

Proposition 16 Assume that the function I(x,y) = h{x) Y\ i=x {y — gi{x)) ai , where 
q ; 6 C and gi{x) are algebraic functions, and h(x) is such that h'(x)/h(x) is a 
rational function, satisfies that 

dl dl 

P{x, y)^-{x, y) + Q{x, y)^-{x, y) = k(x, y) I{x, y), 
ox ay 

where k(x,y) is a polynomial, then I{x,y) is a Darboux function. 

Proof. By Proposition 7 in [H], we have that each gi{x) is a particular solution 
of equation which is algebraic by hypothesis. Therefore, by Theorem we 
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deduce that the irreducible polynomial f(x,y) such that f(x,g(x)) = gives rise 
to an invariant algebraic curve of system (JTJ). For each gi{x) there is an invariant 
algebraic curve f(x,y) = 0, but each invariant algebraic curve f(x,y) = can 
implicitly define several functions gi(x), as much as the degree s of f(x, y) in y. We 
denote by gj t (x) all the algebraic functions defined by the same invariant algebraic 
curve fj(x,y) = 0. Assume that the finite set {gi(x) : i — 1, . . . ,£} defines a total 
of r invariant algebraic curves fj(x,y) = 0, j = l,2,...,r. Hence, the function 
I(x, y) is written as I(x, y) = h(x) Yi(y — 9ji{x)) aji where the product is taken over 
all the possible subindexes. Some of the can be null in this denomination. 

By Proposition |Hl we have that each particular solution y = gj^x) has an associ- 
ated quasipolynomial cofactor, which will be denoted by Mj.(x,y). Each invariant 
algebraic curve fj(x,y) = has a polynomial cofactor denoted by kf.(x,y). We de- 
fine Mo(x, y) := h'(x)P(x, y)/h(x), which is a rational function by hypothesis. Each 
polynomial fj(x,y) decomposes, by TheoremEl as fj(x,y) = hj(x) YliLiiv ~ 9ji( x )) > 
where hj(x) is a function of x such that h'j(x)/hj(x) is a rational function, g^{x) 
are algebraic functions and Sj is the degree of fj(x,y) in y. We define Mj (x,y) : = 
h'j(x)P(x,y)/hj(x) which is a rational function. We have that: 



P(X:V fA +Q{X:y fA 



Pix^h'^x) Yliy-gj^x)) + hj\ 



X 



1=1 



v=l,v^i 



i=l 



M j0 (x, y) + J2 M ^ x > v) h ^ x ) Il> - 9iMl 

K i=l J i=l 

Si \ 

M j0 (x,y) + J2 M ji(^y)\ fjfay). 



i=i 



On the other hand, 



Qf. Qf. 

P{x,y)-^{x,y) +Q(x,y)^- (x,y) = k fj (x, y) fj(x, y), 



dy 



for being fj(x, y) = an invariant algebraic curve with cofactor kfjx, y). Therefore, 
we deduce the following identities: 



M jo (x,y) + ^M jt (x,y) = k fj (x,y), 



i=l 



for j = l,2,...,r. We have that each Mj t (x,y) is a rational function of x, a 
polynomial in gjX x ) an d a polynomial in y, by Proposition 7 in We can change 
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the powers of gj { {x) which are equal to or higher than Sj to a linear combination 
of lower powers by using the expression fj(x, gjX x )) = 0. Therefore, equating 
the powers of y (after changing all the powers gjX x ) to combinations of gj t (x) u , 
< v < Sj), we get that the r identities (fTHjl give a total of r(m — 1) linear 
combinations of the powers gjX x Y \ < v < Sj, where m is the degree of system ((Tj) 
in the variable y. 

We consider the identity 

dl dl 
from which we deduce that 

Sl S2 s r 

Mo + J2 <* U M U + «2,M 2i + + a n M n = k - (19) 

i=l i=l i=l 

As before, we can change each power g^ixY with v > Sj to a linear combination 
of the powers g^ixY \ < u < Sj. And from (fT^j) we deduce an identity which 
is a linear combination of the powers gj^x)", i = 1,2, ... ,Sj, < v < Sj and 
j = l,2,...,r. 

We notice that if aj i = f3j, where j3j G C, for all 1 < i < Sj, and 1 < j < r, 
then the relations given by (fTK|l make the identity (fTTH) compatible with the fact 
that k(x,y) is a polynomial in both variables x and y. Assume that this is not 
the case. Assume that we have otj v ^ aj v for certain v, v. We can assume that 
j = 1 without loss of generality. We can consider the Sj symmetric polynomials 
defined by each fj(x,y) which are linear combinations of the powers gj^x)", i = 
1,2, Sj, < v < Sj, j = 1,2, ... ,r. By using the elimination theory, we can 
eliminate all the appearances of gjX x ) with j > 1 in (|19j) . We obtain in this way 
a relation only involving gi^x). By using the symmetric polynomials associated 
to fi(x,y), we eliminate all the gi^x), except one, which may be gi^x). The 
resulting relation R(x, gi^x)) = can only be of two forms: either R(x,y) is a 
multiple of f\(x,y) or it is not. In the first case, we have that the relation given 
by (j!9)l is a combination of the symmetric polynomials, which are symmetric with 
respect to (71,(2), % = 1, 2,...,s 1 . This symmetry implies that oij v = aj v , for all 
v and v. In the second case, we would get that R(x, y) is a polynomial such that 
R(x, gi^x)) = and the degree of R(x,y) in y is lower than s±. We recall that we 
have already substituted all the appearances of powers of g^ (x) of higher degree by 
the corresponding expression given by the equation f±(x, gi^x)) = 0. The existence 
of a polynomial R(x,y) such that R(x, g^ (x)) = and the degree of R(x,y) in 
y being lower than Si is a contradiction with the fact that j\(x,y) is the only 
irreducible polynomial satisfying fi(x, gi 1 (x)) = 0, modulus associates. 

Hence, we conclude that the only possibility is that = (3j, where (3j G C, for 
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all 1 < i < Sj, and 1 < j < r. We have that: 

I(x,y) = h(x) ^ =Mx)Pi 

and we define h(x) = h(x) YFj=i hj(x)~^. We notice that since h(x) and hj(x), 
j = 1, 2, . . . , r satisfy that its logarithmic derivative is a rational function, h(x) also 
satisfies that h'(x)/h(x) is a rational function. By integration, we obtain that h(x) 
is a Darboux function. We deduce that I(x, y) is equal to h(x) Ylj=i fj( x i v)^ which 
is a Darboux function, as we wanted to show. | 

Example. We are going to describe an example of Proposition so as to 
make the proof clearer. Let us consider the following planar polynomial differential 

system: 

x = — 5 — 5x + 15y 2 — 6x 2 y + lAxy 2 — 9xy 4 , , , 

y = 5 + 2x - 3y - 2xy 2 + 6y 3 - 3y 5 . ^ ' 

This system exhibits two invariant algebraic curves of degree 3: f\(x,y) = with 
fi(x,y) := y 3 — y — x and f 2 (x, y) = with f 2 (x, y) '■= xy 2 — x — 1. Their cofactors 
are, respectively, kf^x.y) = —3(l + 2xy—Ay 2 + 3y 4: ) and kf 2 (x, y) = 5 ki(x, y)/3. We 
factorize the polynomial fi(x,y) as hi(x)(y — gx 1 (x))(y — gi 2 (x))(y — gi 3 (x)) where 
hi(x) := 1 and gi^x), i = 1, 2, 3, are the corresponding y-roots of f\(x, y). It is easy 
to see that the polydromy order of g^ (x) is 1. We perform the same computations for 
f 2 {x,y) and we have that f 2 {x,y) = h 2 (x) (y - g 2l (x) ) (y - g 22 (x) ) where h 2 (x) := x, 
g 2l (x) := y/l + l/x and g 22 (x) := — y/T+ TJx. It is easy to see that the polydromy 
order of g 2i (x), i = 1, 2, is 2. We have that y—g^x) are algebraic particular solutions 
and we can compute their corresponding quasipolynomial cofactors Mj^x^y) which 
are: 

M h (x,y) := ^ T —(-l + 5y + 2y 2 -y 4 + (5 + 2x-2y + 2y 3 )g u (x) + 

+(l + 2zy-4y 2 + 3y 4 K 2 i (x)) ) 

M 2i (x, := — — (by — lOx — 4a; 2 + 6xy — 15y 3 — 6xy 3 + 

2xg 2i (x) 

+ (5 — 15y 2 — Ax 2 y + Qxy 2 — 6xy )g 2i (x)) . 

We have that gi^x), i = 1,2,3 are the y-roots of fi(x,y) and, hence, they satisfy 
the following relationships, given by the symmetric polynomials on the y-roots: 



ti 2 
ti 3 



9ii(x) +9i 2 (x) +9i 3 (x) =0, 
^liC^PiaC^) + 9ii(x)gi 3 (x) + g l2 (x)g l3 (x) = -1 
g ll (x)g l2 (x)g l3 (x) = x. 
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In the same way we have that V 2l : = #21 0*0 + 92 2 ( x ) — and V 22 := g 2l (x)g 22 (x) = 
— 1 — 1/x. These relationships give that M ll +M l2 +M l3 = k\ and M 2o +M 2l +M 22 = 
k 2 where M 2o = h' 2 (x)x/h 2 (x) = (—5 — 5x + 15y 2 — 6x 2 y + lAxy 2 — 9xy A )/x. We 
consider a function of the form: 

i(x, y ) = ^(y-Pi!^)) 01 ^^-^^^) 01 ^^-^^)) 01 ^^-^^)) 1121 ^-^^))^, 

where is such that h'(x)/h(x) is a rational function (we define M (x,y) : = 
h'(x)x/h(x)), aj t G C and gjX x ) are the functions defined above. We also assume 
that X (I(x,y)) = k(x,y)I(x,y) where k(x,y) is a polynomial in x and y of degree 
at most 4. We have that: 

M + o^M^ + a l2 M l2 + a l3 M l3 + a 2l M 2l + a 22 M 22 = k. 

This identity gives rise to five equations relating the g^x) and when equating 
the coefficients of different degrees of y. The equation corresponding to the degree 
4 in y is equal to: 

( h'(x)\ 
— 3 [ cki, + ol\ 2 + ai, + a 2l + a 22 + 3x , . . = fc 4 (x), (21) 
V " Kx)J 

where kj(x) is the coefficient of degree j in y of /c(x, y). We note that /^(x) is a real 
number since k(x, y) is a polynomial in x and y of degree at most 4. Then, equation 
(|2Tj) implies that /i(x) = x b with 6 = — (ai 1 + cti 2 + ai 3 + ct 2l + a2 2 )/3 — k^/9. The 
equation corresponding to the degree 3 in y gives: 

/ 6a u g u (x) \ _ 3(5 + 2x) f a 2l a 2 , 2 \ = ( 
^\l-3q,Ax) 2 2x \q 9 Jx)^q 9 Jx) 31 h 



i=l 



The following equations correspond to the degrees 2, 1 and in y, respectively: 
,h'(x) A (u x .(2g x ix) 2 - 1)\ 3(2x - 5) , s I/N 



^/ 5-2^ 1 .(a;) + 2^ 1 .(a;) 2 \ >A / 



5 + 6x — Ax g 2i (x) 

OL 2 .~ 



2xg 2i (x) 



2 h'(x) 



= 6x - . , + fci(x), 
h[x) 

*f -l + (2x + 5)g u (x)+g u (x) 2 \ , A/ 5g 2i (s) - 10s - 4x 2 
= 5(1 +x)gg + *,(*)■ 
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We consider a common denominator in each one of these equations and by us- 
ing elimination theory among the numerators of these equations and the poly- 
nomials Vj v we deduce that the only possibilities are ai 1 = a>i 2 = ai 3 =: f3\, 
a 2l = a 22 =: (3 2 , h(x) = x& and k(x,y) = -(3/3i + 5/3 2 )(l + 2xy - Ay 2 + 3y 4 ). 
Hence, I(x,y) = fi(x, y) /3l f2(x, y) 132 , which is a Darboux function. 

We know that taking y-roots any Darboux function can be expressed in the form 
I(x,y) = exp {h 2 (x)H r k=1 (y - a k (x))/Hj =1 {y - gj(x))j h x {x)Xf i=l {y - gi {x)) a * 
where a, G C, gi(x), (jj(x) and a,k(x) are algebraic functions, hi(x) is such that 
its logarithmic derivative is a rational function and h 2 (x) is a rational function. 
Next Proposition gives the converse of this assertion, that is, we give the conditions 
that a function of the form I(x, y) must satisfy in order to be a Darboux function. 

Proposition 17 Assume that the function 
I(x, y) = exp 



M*) 5H y a f\l h(x) ]J(y - 9l (x)y 
Uj=i(v ~ 9 Ax)) J f = \ 



where cti G C, gi(x), (jj(x) and ai~(x) are algebraic functions, h\(x) is such that 
its logarithmic derivative is a rational function and h 2 (x) is a rational function, 
satisfies that: 

dl dl 

P{x, y)^-{x, y) + Q{x, y)^-(x, y) = k(x, y) I(x, y), 
ox ay 

where k(x,y) is a polynomial, then I{x,y) is a Darboux function. 

Proof. From Theorem 2 in jSj, we have that each gi(x) is a particular solution 
of (0), so it has an associated quasipolynomial cofactor (see Proposition [SJ) which 
we denote by Mi(x,y). We define $(x, y) := exp{h 2 (x) Ax(x, y)/A (x, y)} where 

Ai(x,y) = YTk=i(y ~ a k(x)) and A (x,y) = nj=i(z/ ~ 9j( x ))- Since <*C0 = kI , we 
deduce that 

x ($(x, y )) = (k( x , y ) - j^j^fo y)-Jl a * M ^ y)\ ^ y)- 

Therefore, we are under the hypothesis of Proposition El and Theorem El Hence, 
we realize that the proof of this assertion goes exactly as the proof of Proposition 

m i 

As a consequence of Theorem El and Propositions EJ El El an d El we can 
establish the following result. 
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Theorem 18 Assume that system (JTJ) has a first integral or an integrating factor 
of the form 



where ai G C, Qi{x), cjj(x) and at{x) are algebraic functions, h(x) and h\(x) have 
a rational logarithmic derivative and h2(x) is a rational function. Then, I(x,y) is 
a Darboux function. 

Proof. We are under the hypothesis of Proposition [Tol or IT71 In this case, k(x,y) is 
identically zero or k(x,y) is minus the divergence of the system. We deduce that 
I(x, y) must be a Darboux function. | 
We notice that when we apply the method of constructing first integrals given in 
[H] and all the h(x) and g%{x) are completely determined functions, by Theorem fTSl 
we have that this first integral is a Darbouxian function. However, the reciprocal is 
not true. As some examples in jB] show, we may have a system with a Darbouxian 
first integral, but when we apply the method described in p] we get a nonlinear 
superposition principle. That is, not all the gi(x) introduced in the ansatz I(x, y) are 
determined. If we choose those undetermined Qi{x) as algebraic particular solutions, 
we will have that the first integral given by the superposition principle becomes a 
Darboux function. 
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